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We study different ways of symmetry fractionalization in Z 2 spin liquids on the triangular lattice. 
Our classification can be used to identify the symmetry fractionalization in the Z 2 spin liquid 
reported in recent density-matrix-renormalization-group simulations for J 1 -J 2 spin model on the 
triangular lattice. We find 64 types of symmetry enriched Z 2 spin liquid states on triangular lattice. 
Besides 8 states constructed in Schwinger-boson parton wavefunctions, 12 more states can be realized 
in Abrikosov-fermion parton construction. Within a larger gauge group than SU(2), the rest 40 
states are also found in a spin-3/2 system. Among 20 types of Abrikosov-fermion Gutzwiller- 
projected wavefunctions, No.B5 state is a promising candidate for the Z 2 spin liquid for J 1 -J 2 spin 
model on the triangular lattice. No.BS lies close to Dirac spin liquid (DSL). However, variational 
Monte Carlo simulation find that DSL has a good variational energy and J 1 -J 2 spin model cannot 
open a gap for spinons on top of DSL to stabilize No.B5 state. 


I. INTRODUCTION 

Together with global symmetries, Long-range entan¬ 
glement has been playing an important role in charac¬ 
terizing and classifying the many-body wavefunctions. 
While the former describes various symmetry-breaking 
phased, the long-range entanglement precisely rephrases 
the intrinsic topological order in the gapped quantum 
phases of mattePEI, as it always survives under local 
unitary operations without closing the gajP. This topo¬ 
logical order is fully characterized by the quasi-particle 
(anyon in the bulk) statistical^ (and also the chir al cen - 
tral charge if there are chiral gapless edge stateif^^) . 
In two spatial dimensional systems, the anyon statistics 
and edge chiral central charge (up to E$ bosonic quan¬ 
tum Hall states) is fully described by the modular S and 
T matrices that generate a projective representation of 
SL( 2, Z) on the degenerate-ground-state Hilbert space on 
a torusPnnSEIIl 

When both of the global symmetries and long-range 
entanglements are present, symmetries further enrich the 
classifications of the gapped many-body wavefunctions 
with intrinsic topological orders such that anyon excita¬ 
tions in topological ordered states can carry fractional¬ 
ized quantum numbers of symmetries 20 22 . For exam¬ 
ple, v = 1/3 Laughlin state has anyons carrying the 
electric charge e/3 which is the fractionalized quantum 
number of electromagnetic 17(1) symmetry.® The quan¬ 
tum spin liquid is another concrete example 23 ^. In this 
work we are primarily interested in spin liquids with the 
Z 2 topological order. Besides the trivial particle excita¬ 
tion, Z 2 spin liquid supports three types of anyon ex¬ 
citations: the bosonic spinon b, the vison v and the 
fermionic spinon /, which is the bound state of b and 
v. Spi nons in Z 2 spin liquid state carry a half-integer 
spii J — I—I representing the fractionalization of SO(3) spin 
rotation symmetry. Moreover, anyons (6, v and /) can 
also carry fr action alized quantum numbers of the crys¬ 
tal symmetry® 0 21 \ e.g., the anyons can carry half-integer 


angular momentum or the translations in two different 
directions can become anticommutative when acting on 
an anyorP3£ 2 7 Different ways of symmetry fractionaliza¬ 
tion distinguish different types of Z 2 spin liquid states, 
and this difference can be detected by experiments and 
numerical simulation ;] 2 d 25 ! ^®. 

The search for quantum spin liquids is a long sought 
goal in condensed matter physic^Ml, Z 2 spin liquid 
states have been discovered by numerical studies in sev¬ 
eral two-dimensional frustrated spin systemiP®®^, anc [ 
recent ly in the Ji-J 2 Heisenberg model on the triangular 
lattic^MSJ. To explore the symmetry enriched proper¬ 
ties, we classify different ways of fractionalizing crystal 
and time-reversal symmetries on the triangular lattice 
in Z 2 spin liquid state. Taking into account two physical 
constraints that the spinons (b and /) carry a half-integer 
spin and the vison is described by an odd Z 2 gauge the- 
0 r yl 36 l' 37 l, we find that there are in total 64 different ways 
to fractionalize the crystal and time-reversal symmetries, 
each of which describes one type of Z 2 spin liquid state. 
These 64 states can labeled by six Z 2 quantum num¬ 
bers labeling commutation relation fractionalization and 
quantum number fractionalizatiorP^, and we propose fu¬ 
ture numerical studies to measure these quantum num¬ 
bers with methods described in Ref. 12711271 

The wavefunctions of Z 2 spin liquids can be written 
in terms of spinon (6 or /) representations. Actually, 
a physical spin is fractionalized into either two bosonic 
spinons ( b ) or two fermionic spinons (/), kno wn as the 
Schwinger bosorP^ and the Abrikosov fermiorpS1221, re¬ 
spectively. We can write down the mean held wavefunc¬ 
tions for bosonic (6) and fermionic (/) spinons and im¬ 
plement Gutzwiller projection to project out double oc¬ 
cupied configurations in the mean held wavefunctions. 
The symmetry fractionalization can be studied using 
the projective symmetry group (PSG) analysis on the 
Gutz willer- projected wavefunctions (GPWF) for Z 2 spin 
liquid^ 0 21 . Out of these 64 states, eight has been con¬ 
structed using Schwinger-boson parton construction in 
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Ref. 32| and EH We find that using Abrikosov-fermion 
parton construction with an SU(2) gauge symmetry, 
12 more states can be constructed besides these eight 
states. The connection between Schwinger-boson and 
Abrikosov-fermion constructions is established for the 
eight GPWFs. For Abrikosov-fermion GPWFs, varia¬ 
tional Monte Carlo method is implemented to search for 
the candidate Z 2 spin liquid for J 1 -J 2 spin model on the 
triangle lattice. No.B5 state in TABLE [TT] that lies in the 
neighbor of U(l) Dirac spin liquid (DSL) is promising. 
However, our preliminary variational Monte Carlo sim¬ 
ulations find that the spinon gap on top of DSL is not 
favored and DSL has a good variational energy for the Ji- 
J 2 spin model, consistent with previous systematic vari¬ 
ational Monte Carlo seardPH where gapless spin liquid 
is reported. The parton construction only with a SU(2) 
gauge symmetry doesn’t give rise to all possible ways of 
symmetry fractionalization in Z 2 spin liquid^ 2 -’ 22 . We 
consider a spin-3/2 system and find out the rest of the 
44 states within a gauge symmetry larger than SU(2). 

The rest of the paper is organized in the following: In 
Sec. [II] we use projective method to investigate Z 2 spin 
liquids, including a general consideration ahead, then the 
classification of 20 kinds of fermionic spinons on triangu¬ 
lar lattice and mapping between Schwinger boson states 
to Abrikosov fermion states bridged by vison’s fraction¬ 
alized quantum numbers derived from a dual Z 2 gauge 
held theory. In the next Sec. Ill we carry out prelim¬ 
inary variational Monte Carlo simulations on some of 
these states. Enlightened by the discussion in Sec. |TT| 
in Sec. [IVjwe propose to construct all the Z 2 spin liq¬ 
uid states beyond conventional Abrikosov fermion and 
Schwinger boson formulations. 


II. SYMMETRY FRACTIONALIZATIONS IN Z 2 
SPIN LIQUIDS 

A. General Consideration 

We begin by classifying different ways to fractional¬ 
ize crystal and time-reversal symmetries for one type of 
anyon. In Z 2 spin liquid, an anyon can carry fractional¬ 
ized quantum numbers. This is because a physical wave 
function always contains pairs of anyons (or in general 
anyons that fuse into a trivial excitation). In presence of 
symmetries, the wavefunction carries a linear representa¬ 
tion of the symmetry group G. Due to fractionalization, 
anyons carry a projective representation of G, which fea¬ 
tures fractionalized quantum numbers of the global sym¬ 
metry. Precisely, the relation f(g)f(h) = f(gh ) for a 
linear representation / acting on two group element g 
and h can be fractionalized into f(g)f(h ) = to(g, h)f(gh) 
in a projective representation, where the factor uj(g , h) 
labels the symmetry fractionalization. 

The symmetry fractionalization embedded in the pro¬ 
jective representation must be consistent with the anyon 
fusion rule. Particularly since two anyons of the same 


type in Z 2 spin liquid fuse into a trivial excitation, the 
factor uj(g,h) can only be ±1. Furthermore, two projec¬ 
tive representations that differ by a linear representation 
represent the same topological order, because we can al¬ 
ways attach a trivial local excitation carrying a linear 
representation to an anyon. Therefore if we consider only 
one type of anyon, different symmetry enriched topolog¬ 
ical orders are classified by inequivalent projective rep¬ 
resentations with co(g,h) = ±1, and mathematically this 
is given by the second group cohomology with Z 2 coeffi¬ 
cient, h 2 (g,z 2 )P 

Different elements in FI 2 (G, Z 2 ) are labeled by a set 
of uj factors which indicates how an algebraic relation 
gh - ■ ■ = 1 in the group algebra of G is promoted 
to f(g)f(h) ■ ■ ■ = ±1 in a projective representation. 
These w factors are the cocycle variables in the lan¬ 
guage of group cohomology. Physically they indicate 
how the action of symmetry group is fractionalized on 
the anyon. For example, two commuting symmetry op¬ 
erations XY = YX can become anticommuting XY = 
—YX when acting on an anyon, and a relation X n = 1 
can become X n = — 1 when acting on an anyon, indi¬ 
cating that the anyon carries a fractionalized quantum 
number of ,V 

On the triangular lattice (Fig. [Tb, the symmetry group 
we consider here is G = G p Q m x Zj, where G P Q m and Z 2 
denotes the space group of the lattice and time-reversal, 
respectively. For this symmetry group, the second group 
cohomology can be calculated using the HAP package® 
in the GAP systeuP^. It can be shown that when cal¬ 
culating the group cohomology, the subgroup of transla¬ 
tions can be separated from the point group Cq v and the 
time-reversal symmetry, 

% 2 (G,Z 2 ) = Z 2 x 'H 2 {C(, v x Z^,Z 2 ), (1) 

where the first Z 2 denotes the fractionalization of the 
commutation relation between the two translation oper¬ 
ations, T 1 T 2 = ±T 2 Ti. Hence when classifying projective 
representations on the triangular lattice, the translation 
symmetry can be considered separately from the point 
group and time reversal symmetries 4 -. (This separation 
is a special property of the triangular lattice space group 
and may not apply to other lattices, e.g., the square lat¬ 
tice). 



FIG. 1. Triangular lattice in a oblique coordinate system. 

Furthermore, the symmetry fractionalization of the 
point group and time-reversal symmetries is classified by 
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six more Z 2 labels as indicated by 'H 2 {Cq v xZ 2 , Z 2 ) = Z 2 . 
In this paper we choose the following six Z 2 variables 
Wjj, ujj, lot, uicrti Wpt, where each uj (X = cr, g,I) labels 
a quantum number fractionalizatiorP^ X 2 = ±1 for the 
symmetry operation X , as listed in Table [I] Here a and g 
are two mirror reflection symmetry operations as labeled 
in Fig. [T] T is the time-reversal symmetry and / is the 
inversion (or 180-degree rotation) symmetry / = (ag) 3 . 
Together with the translations Tf and T 2 they generate 
the space group G p em and their action is the following in 
the coordinate system shown in Fig. [TJ 

T\ : (ri,r 2 ) ->■ (rq + l,r 2 ), 

T 2 ■ (ri,r 2 ) ->■ (ri,r 2 + 1), 

M : (ri,r 2 ) -> (rq - r 2 ,-r 2 ), (2) 

o : (ri,r 2 ) ->■ (ri,ri - r 2 ), 

I : (ri,r 2 ) ->■ (-rq,-r 2 ). 

The symmetry fractionalization of one type of anyon 
is classified by the aforementioned seven Z 2 labels, as 
T~L 2 (G, Z 2 ) = Z 2 . In general each type of anyon in the Z 2 
topological order can have different symmetry fraction¬ 
alization, but their symmetry fractionalization must be 
consistent with the anyon fusion rule, so the symmetry 
fractionalization of two anyons determines the symmetry 
fractionalization of the third. Moreover, some combi¬ 
nations have obstructions and can only be realized on 
the surface of a three-dimensional nontrivial symmetry- 
protected topologically ordered state. Here we further 
restrict ourselves to Z 2 spin liquid states that satisfies 
the following two physical conditions. First, we assume 
that the vison is described by an odd Z 2 gauge theory 
on the dual lattice and therefore has a fixed way of sym¬ 
metry fractionalization, which is studied in Sec. |II C| As 
a result, there is a one-to-one correspondence between 
the symmetry fractionalization of bosonic and fermionic 
spinon. Second, we assume that the spinons carry a half¬ 
integer spin and therefore T 2 = —1. This fixes the label 
lot = — 1 for the bosonic and fermionic spinons. With 
these two constraints, the symmetry fractionalization in 
a Z 2 spin liquid state is labeled by six Z 2 labels, and 
therefore there are 2 6 = 64 different types of spin liquids. 

The symmetry fractionalization quantum numbers we 
list in Table [I] can be directly measured using the sym¬ 
metry quantum numbers and symmetry representation of 
Schmidt eigenstates in the entanglement spectrum of the 
ground state wave functions obtained in numerical simu¬ 
lations including the DMRG studies 25 2 ‘. First, the com¬ 
mutation relation fractionalization number u 12 can be 
measured from the ground state momentum in the direc¬ 
tion of T 2 on a torus with odd number of unit cells in the 
direction of T\. Second, the fractional quantum number 
of the crystal symmetries u) p and ujj can be measured 
from the symmetry quantum numbers of the ground state 
wave functions. Lastly, the fractional quantum numbers 
X 2 = ojx can be measured as a one-dimensional SPT 
invariant by calculating how the symmetry acts on the 
Schmidt eigenstates in the entanglement spectrum, as all 


TABLE I. Quantum numbers labeling different projective rep¬ 
resentations of the symmetry group H 2 (G, Z 2 ). 


No. 

Algebraic relation 

Quantum number 

1 

T x T 2 = ±T 2 Ti 

0 J 12 

2 

g 2 = ±1 


3 

a 2 = ±1 

UJ(j 

4 

1 2 = ±1 

OJI 

5 

T 2 =±1 

lot 

6 

(gT) 2 = ±1 

COfit 

7 

( 0 -T ) 2 = ±1 

bJcrt 


symmetry operations listed in Table [T] can be realized ei¬ 
ther as an on-site antiunitary operation or an unitary and 
left-right-exchanging operation if we view the system as 
a quasi-one-dimensional system. 


B. Symmetry fractionalization in parton 
constructions of Spin Liquids 


Z 2 spin liquid can be described by GPWFs in parton 
construction approaches. In these approaches the spin 
operator is decomposed into two spinons 


Si = 


1 

2 


0’i a (T apO’iPi 

otfi 


O'icx. 


ficn bia 


( 3 ) 


where a ia can be a fermionic spinon f ia or a bosonic 
spinon 6j a in the Abrikosov fermion and Schwinger bo¬ 
son constructions, respectively. Spinons have the one- 
particle-per-site constraint Y^ a a ia aia = 1- Both parton 
constructions have an emergent gauge symmetry: In the 
Schwinger boson construction the physical spin operator 
in Eq. © is invariant under a U(l) symmetry transfor¬ 
mation bi a —► e lSi bi ai whereas in the Abrikosov fermion 
construction it is invariant under an SU(2) gauge trans¬ 
formation ifji gipi, where ip = and g is an 

element of the SU(2) group generated by Pauli matrices 

rW 3 . 

Since we assume that the spinons carry a half-integer 
spin, T 2 = — 1 is fixed. Hence the spinon transforms 
under T as Tai a = e a pa,ip. For other symmetry group 
actions, in the PSG language the symmetry action X 
is promoted to GxX, where Gx is an element of the 
gauge group (SU(2) for Abrikosov fermions and U(l) for 
Schwinger bosons). Using the generators of the space 
group, the generic PSG solution for spinons is given as 

gI (g, r 2 ) = wg, Gj,* (n, r 2 ) = 1, 

G^(r 1 ,r 2 )=^ +r2(r2+1)/ V, W 

G b /(n,r 2 ) =w ^2 ( ’ 1+1)/2 .9,t, 


where g^^a are independent of the site coordinates and 
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satisfy the algebraic relations 

9l = gl = w CT , (3 m3<t) 6 = wj, ^ 

- -ffcrfiC = w CT r- 

For fermionic spinons, g^,g a £ SU(2) while g^.,g a are 
17(1) phases for bosonic spinons. Different Z 2 spin liquids 
can be constructed by solving these equations with differ¬ 
ent combinations of symmetry fractionalization quantum 
numbers. For Abrikosov fermion constructions with an 
SU(2) gauge symmetry, we find that 20 different spin liq¬ 
uid states can be realized, as listed in Table |ll| We note 
that on the kagome lattice which has the same space 
group, Abrikosov fermion construction also realizes 20 
different stateiP^, and more states can be realized on 
an anisotropic triangular lattic^®, which has a smaller 
space group. On the other hand Wang and VishwanathES 
showed that the Schwinger boson construction with U(l) 
gauge symmetry can realize eight Z 2 spin liquids states. 

We note that neither Schwinger boson nor Abrikosov 
fermion constructions can realize all 64 different com¬ 
binations of symmetry fractionalization. This is a re¬ 
striction due to the gauge symmetry presented in these 
parton constructions. As an example, in the Abrikosov 
fermion construction if = —1 and = +1, u>i must 
equal to — 1 because g a = r° without any other choice 
for the algebra spanned by {r°, ir 1 , ir 2 , ir 3 } in regard to 
the two dimensional representation of SU(2) group. The 
Schwinger boson construction is more restrictive because 
the U(l) gauge group is smaller than SU(2). The other 
44 states that are not constructed in this section requires 
a parton construction with a gauge symmetry larger than 
SU(2) and we will give an example of their construction 
in Sec. ITVl 

The PSG analysis discussed in this section natu¬ 
rally describes the symmetry fractionalization of bosonic 
spinons and fermionic spinons in the Schwinger boson 
and Abrikosov fermion constructions, respectively. How¬ 
ever, the symmetry fractionalization of the other two 
types of anyons are not revealed explicitly in this anal¬ 
ysis. In the next section we shall see that in all these 
parton constructions the vison has the same symme¬ 
try fractionalization quantum numbers and the symme¬ 
try fractionalization of the other type of spinon can be 
derived by combining the quantum numbers of spinon 
and vison. Therefore there is a one-to-one correspon¬ 
dence between bosonic and fermionic PSG solutions, and 
we see that the 20 different spin liquid states we con¬ 
struct using Abrikosov fermion covers the eight states 
constructed with Schwinger boson. The connection be¬ 
tween fermionic and bosonic states is listed in Table ITTT1 


C. Vison as a bridge between fermionic and 
bosonic spinons 

We start in this subsection to investigate the fractional¬ 
ization of vison on triangular lattice using a gauge theory 


TABLE II. All twenty kinds of Abrikosov fermions on trian¬ 
gular lattice. For states No. A1 - No. A10, W 12 = 1; for No. 


B1 - 

No. 

BIO, 

Cdl2 

= - 

-1. 







No. 


Ld(j 

Wl 

Uflt 

Id at 

Se 

g<y 

U a 

up 

Ury 

a 0 

A1 

+ 

+ 

+ 

— 

— 

T° 


T 1 ’ 3 

T 1 ’ 3 

T 1 ’ 3 

r 1 ’ 3 

A2 

+ 

— 

— 

— 

+ 

T° 

ir 2 

0 

0 

0 

0 

A3 

+ 

— 

— 

— 

— 

T° 

ir 1 

T 1 

T 1 

T 1 

T 1 

A4 

— 

+ 

— 

+ 

— 

ir 2 

T° 

0 

0 

0 

0 

A5 

— 

+ 

— 

— 

— 

ir 1 

7~° 

T 1 

T 1 

T 1 

T 1 

A6 

— 

— 

+ 

+ 

+ 

ir 2 

ir 2 

0 

0 

0 

0 

A7 

— 

— 

+ 

— 

— 

ir 1 

ir 1 

T 1 

T 1 

T 1 

T 1 

A8 

— 

— 

— 

+ 

— 

ir 2 

ir 1 

0 

0 

0 

0 

A9 

— 

— 

— 

— 

+ 

ir 1 

ir 2 

0 

0 

0 

0 

A10 

— 

— 

— 

— 

— 

ir 1 

ir 3 

0 

0 

0 

0 

B1 

+ 

+ 

+ 

— 

— 

T° 

T° 

0 

0 

T 1 ’ 3 

T 1 ’ 3 

B2 

+ 

— 

— 

— 

+ 

T° 

ir 2 

0 

T 1 ’ 3 

0 

0 

B3 

+ 

— 

— 

— 

— 

T° 

ir 1 

0 

T 3 

T 1 

T 1 

B4 

— 

+ 

— 

+ 

— 

ir 2 

T° 

T 1 ’ 3 

0 

0 

0 

B5 

— 

+ 

— 

— 

— 

ir 1 

7~° 

T 3 

0 

T 1 

T 1 

B6 

— 

— 

+ 

+ 

+ 

ir 2 

ir 2 

0 

0 

0 

0 

B7 

— 

— 

+ 

— 

— 

ir 1 

ir 1 

0 

0 

T 1 

T 1 

B8 

— 

— 

— 

+ 

— 

ir 2 

ir 1 

T 1 

0 

0 

0 

B9 

— 

— 

— 

— 

+ 

ir 1 

ir 2 

0 

T 1 

0 

0 

BIO 

- 

- 

- 

- 

- 

ir 1 

ir 3 

T 3 

T 1 

0 

0 


and PSG analysis. As a matter of fact, vison is a little 
bit different from fermionic spinon and bosonic spinon. 
There is only one kind of vison here thus it is meaning¬ 
less to talk about classification. oj ] j 2 = — 1 because a vison 
sees a spin-1/2 quanta. As mentioned in the introduction 
section, vison is nothing but the vortex excitation in an 
odd Z 2 gauge theory 36 37 . To begin with a common Z 2 
gauge theory, 

ff = ^E^- r EL(6) 

* *(») 

where variables t x = ±1 are placed on bonds of trian¬ 
gular lattice and i runs on the dual honeycomb lattice. 
The conserved Ising gauge charge n 1(a) T~l(a) ^ ^ a7 * 
odd Ising gauge theory means the product measurement 
of all six lines emanating from site a, which also can be 
understood in the quantum dimer model where there is 
only one dimer on each sit^S. By defining new opera¬ 
tors on the dual honeycomb lattice erf = T i(i) an d 
of = Ilz<i T i a F u ^y Frustrated Ising Model (FFIM) on 
dual lattice can be written as®®] as 

T^dual = 7 d ^ ' ^ij&i ®j b p ^ ^ C if, (7) 

(*J> * 

in which there is exactly one = —1 along each 
hexagon, thereby the fractionalized quantum numbers for 
vison may be calculated from this FFIM. In the dual Ising 
language, the Ising spins cr z are vison creation operators. 

As shown in Fig. [2j the honeycomb lattice has the 
same space group symmetry as the triangular lattice, 
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(a) 


FIG. 2. Fully frustrated Ising model on dual honeycomb lat¬ 
tice of triangular lattice. Red and bold bonds represent frus¬ 
trated ones, namely negative while otherwise are all positive. 
A and B denote two sublattice for the dual honeycomb lattice. 


which acts in the following way 


Ti : (n, r 2 , a) (n + 1 ,r 2 ,a), 

T 2 : (rq, r 2 , a) -)• (ri,r 2 + l,ct), 

. . , 1 + 3a 1 . 

M : (?T,r 2 ,a) -A (rq - r 2 --—, -r 2 - -,-a), ^ 

a : (ri,r 2 , a) (r i,ri - r 2 ,a), 

2 1 

I : (rq,r 2 ,a) ->• (-n - ^,-r 2 - -,-a). 

where a = ±1 denotes the two sublattices A and B in 
a unit cell (see Fig. [2]). Note that certain symmetry 
operations exchange the sublattice index in Eq. [8] To 
be simple and convenient, we define rq = rq + IT and 
r 2 = r 2 — | to represent integer coordinates for hexagons 
in a same coordinate with triangular lattice. Here we 
derive a specific PSG solution for vison to compute all 
the corresponding gauge invariant fractionalized quan¬ 
tum numbers. It is learned that a U(l) gauge symmetry 
as used in Ref. EH is enough to express the state as il¬ 
lustrated in Fig. [2j After a bit calculation, we get the 
explicit PSG solution to vison 

G v Tl (r r,r 2 ,a) = 1, G£ 2 (rq, r 2 , a) = (-l)*i(i-a)( 9a ) 

G” (rq, r 2 , a) = (-) ? i + 5^(f 2 +3) i (i-a)(i+f 2 ) ; ( 9b ) 

G£(n, r 2 , a) = (_)^ + ^TGi+Di(i+«Fi, ( 9c ) 

G?(r 1 ,r 2j a) = (-r i i 1 +“. (9d) 

These PSG solutions can then be used to compute the 
symmetry fractionalization quantum numbers for the vi¬ 
son, and the result is summarized in Table [TTT1 

These results of vison’s symmetry fractionalization 
complements the classification of spinon symmetry frac¬ 
tionalization discussed in the previous section using par- 
ton constructions, which only give the symmetry frac¬ 
tionalization of one type of anyon (the Schwinger boson 
and Abrikosov fermion construction gives the symmetry 
fractionalization of bosonic spinon and fermionic spinon, 


TABLE III. Mapping rules between fermionic states and 
bosonic states through their fractionalized quantum numbers 
for Z 2 spin liquids on triangular lattice. 


fermion 

vison 

boson 

U 12 

^12 = — 

Ul2 

4 

u V p — + 

4 

4 

U& = + 

ul 

ul 

UJi=- 

u b i 

^ fit 

Upt = + 

u b p,t 

u <Jt 

Uat = + 

ult 


mapping rules 
w 12 — U 12 ' w 12 

4 = -u>l ■ UJ ” 
u\ = -Cjl ■ U V a 

Uj) = —02/ ■ Ul} 
b f v 

Upt = UJit ' 
b f v 

Ucrt = U J at • U at 


TABLE IV. Mapping from bosonic states to fermionic states 
for Z 2 spin liquids on triangular lattice. 


No. 

bosonic states as labeled 
in Ref. ETJwith (pi,p 2 ,P 3 ) 

fermionic states 
as in Table [TT| 

1 

(0,0,0) 

B7 

2 

(1,0,0) 

A7 

3 

(0,1,0) 

B1 

4 

(1,1.0) 

A1 

5 

(0,0,1) 

B5 

6 

(1,0,1) 

A5 

7 

(0,1,1) 

B3 

8 

(1,1,1) 

A3 


respectively). Combining these two results, the symme¬ 
try fractionalization of the third type of anyon can be 
determined. For any quantum number ujx , the values 
of the three anyons are related as oj x = sx<u x uj x be¬ 
cause / is a bound state of b and v, and sx = ±1 is 
a twist factor that arises for symmetry operations that 
exchanges anyon s as a result of the mutual statistics 
between b and / ) 22 l 25 l 26 l. Furthermore, the symmetry 
fractionalization of vison also provides a mapping be¬ 
tween symmetry fractionalization quantum numbers of 
two types of spinons, as summarized in Table |III| The 
already well known eight kinds of Z 2 spin liquids within 
bosonic construction 41 can be mapped exactly to part 
of the states we have found as shown in Table [TTJ where 
the fractionalized quantum numbers for bosonic states 
can be obtained explicitly by means of the PSG solu¬ 
tions given in Ref. ED Here we summarize these results 
in Table |IV[ Also note that this mapping can be inde¬ 
pendently checked using the fact that a nearest-neighbor 

only bosonic RVB wave function Pq exp ^ije^b^b^ |0) 
can be exactly mapped to a Gutzwiller projected BCS 
wave function P G exp Qje a pflfl |0)l51152lif |^.| = 


and <f>p = <pj, + 7r, which means the fermionic plaquette 
flux <j)f differs n from bosonic plaquette flux. 
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III. VARIATIONAL MONTE CARLO RESULTS 

In this section, we will implement VMC simulations to 
search for the candidate Z 2 spin liquid for the J 1 -J 2 spin 
Heisenberg model on the triangular lattice 

H = Ji ^2 Si ■ Sj + J 2 Y, Si ■ S ( 10 ) 

in) m) 

where (ij) and (( ij )) are for the nearest neighbor (nn) 
and 2nd nn bonds, respectively. Previously the VMC 
method has been used to study other spin models on the 
triangular lattice to search for spin liquidpHMI 

Here we 

choose to st udy the J 1 -J 2 as motivated by recent DMRG 
studietPHI. 

Our variational wave function is obtained from the fol¬ 
lowing mean field Hamiltonian for fermionic spinons 

hmf = - y; ip] Ujjtpj+^ M a v , jT°V’j (n) 

ij i 

with the double ip = (/if, /^) T and the mean field pa¬ 
rameters Uij = (ipiipj). 

Let’s consider mean field parameters up to 3rd nn 
bonds (u a , up and m 7 for nn, 2nn and 3nn). The allowed 
non-zero mean field parameters are listed in TABLE [TTj 
Among these 20 different kinds of spin liquids, only No. 
Al, Bl, B3, B5 can be gapped if only mean field param¬ 
eters up to 3rd nn bonds are considered. Because of a 
vanishing mean field parameter on nn bonds, No. Bl and 
B3 are unlikely to be the ground states of J 1 -J 2 model 
on the triangular lattice. No. Al lies closed to U(l) spin 
liquid with large spinon fermion surface while No. B5 is 
in the neighbor of DSL. Without any further parameters, 
DSL is already a good approximation of t he gro und state 
of Ji spin model on the triangular lattice®^. Gapped 
Z 2 spin liquid is most likely to lie in the neighbor of DSL 
and No.B5 is a promising candidate for the Z 2 spin liquid 
on the triangular lattice. We will study the different per¬ 
turbations on top of DSL. Furthermore, the B5 state also 
corresponds to the zero-flux Schwinger-boson mean field 
state studied in Refs. m and HR Explicitly, we consider 
the mean field Hamiltonian 

#dsl = -^sign(y)i/Jr 3 V>j + A* ^2 V’] T Vi 

(ij) i 

+A ^2 sign (ij)ip\r 

(((ij))) i 



FIG. 3. The variational-parameter dependence of local spin 
correlations on nearest neighbor(nn), 2nd nn (2nn) and 3rd nn 
(3nn) bonds. DSL is taken as a reference state and 3 A (Si • Sj) 
is relative local spin correlations. The variational parameters 
are AF Neel spin order h in (a), on-site paring p in (b) and 
3nn bond hopping amplitude A in (c). 


where the signs for mean field parameters on nn and 3nn 

bonds are determined by PSG solutions in Eq. = taken as a reference state. Local spin correlation varia- 

(cos(^(n + r 2 ),sin(^(n + r 2 )))) is the coplanar AF tions are studied on different bonds with increasing h, n 

spin order varying in the 3x3 enlarged unit cell. an< ^ ^ see possible instabilities on top of DSL. 

The VMC is carried out on a 12x12 lattice. DSL (/i = In Fig. [ 3 ] (a), with increasing h, the spin correlations 
A = h = 0) has the local spin correlations 3(S; • Sj) = on nn bonds decreases implying that the AF Neel spin 

—0.5285(2) on nn bonds, 3(S.j • Sj) = 0.2135(9) on 2nn order is favored by J\ spin model on the triangular lat- 

bond and 3(S; • Sj) = 0.0044(8) on 3nn bonds. DSL is tice. When h ~ 0.52, the best variational energy for J\ 
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model is En/{N J{) = —0.54142(7) I 57 * 58 ] Meanwhile, the 
increment of spin correlation on 2nd bonds suggests that 
J 2 term will kill the spin order in agreement with DMRG 
calculation d 34 * 35 ! 

The on-site chemical potential /z gaps the Dirac nodes 
of DSL. With increasing /z, the spin correlations on nn 
and 2nn bonds increases, but very little as shown in Fig. [3] 
(b). The bond parameter A on 3nn bonds is also able to 
open the spinon gap. As shown in Fig. [3] (c), when A is 
small, the spin correlations on nn and 2nd bonds increase 
with increasing A. Therefore, neither /z or A is favored by 
Ji-J -2 spin model on the triangular lattice. The DSL has 
a good variational energy E/{NJ\) = —0.5072(3) when 
J 2 /J 1 = 0.1. From these simple attempts, we find that it 
is not easy to open the spinon gap on top of DSL for J 1 -J 2 
model on the triangular lattice. This is consistent with 
previous VMC simulations in which gapless spin liquid 
is reported^. In addition, we note that similar gapless 
U(l) spin liquid states are also obtained in VMC studies 
on the kagome latticd 521 ^, in contrast to th e gap ped Z 2 
spin liquid state obtained by DMRG studied®^. 


IV. REALIZATION OF Z 2 SPIN LIQUIDS 
BEYOND SCHWINGER BOSON AND 
ABRIKOSOV FERMION CONSTRUCTIONS 


As we have demonstrated previously with a concrete 
example, the Abrikosov fermion construction can not re¬ 
alize all 64 kinds of spin liquid states on triangular lattice 
due to that the parton construction has only an SU(2) 
gauge symmetry and the projective representation is re¬ 
stricted in the linear space of an SU(2) fundamental rep¬ 
resentation. Naturally in this way, this restriction can 
be lifted and more Z 2 spin liquid states are realizable 
if we consider a larger gauge symmetry in the parton 
construction, for example, tensor product of two SU(2) 
fundamental representations. We again use our previ¬ 
ously example = — 1 ,w„ = +1 as an illustration. If 
the gauge group is enlarged to contain products of two 
SU(2) matrices, g a can be chosen as g a = (It 1 ) x (it 1 ) 
and g M = it 3 , therefore resulting in w/ = +1 which is a 
new state beyond Tabic [TT] 

Here we propose a parton construction to realize this 
kind of exotic states. Motivated by the work of Wang and 
XiM we consider a spin-1 system with three degenerate 
orbits, each occupied by one spin-4 electron, and we as¬ 
sume that the spins on three orbits are aligned due to a 
large Hund’s coupling. The spins on the three orbits can 
be fractionalized into three sets of Schwinger fermions, 
s ia = fjaa^aefiad- where a = 1,2,3 denotes the orbital 
degree of freedom. On the other hand, note also that an 
even number of half spin within each unit cell will change 
the fractionalized quantum numbers of vison. For illus¬ 
tration, ujy 2 will change its sign since vison travels around 
a unit cell but does not see net Z 2 gauge charge. 


Using the Abrikosov fermions, the total spin operator 


5 , = Si! + Si 2 + Sis = ^ pi 


iacx. ** 


fifiafi 


a—1 


(13) 


and the orbital occupation number 

n a = Yl fiaafiaoL (14) 

a 

are both invariant under the gauge transformation 
SU(2)xSU(2)xSU(2), where the three SU(2) subgroups 
act on the three copies of Abrikosov fermion opera¬ 
tors, respectively. In other words, if we group the 
Abrikosov fermion operators as a six-component vector 

i’i = (/sit fin /»t fl 2| /at fin) ’ then the § au s e 

transformation acts as the direct sum of three SU(2) ma¬ 
trices, g = gi © g 2 © 33 - 

Next we consider a Abrikosov fermion mean field 
wave function that breaks the gauge symmetry down to 
U(l) xU(l) xZ 2 , where the three subgroups are generated 
by (ir 3 )©(-ir 3 )®l, l©(ir 3 )©(-ir 3 ), and (-1)©1®1, re¬ 
spectively. Although the invariant gauge group (IGGj— 
appears to be larger, this still realizes a Z 2 spin liquid 
because a compact U(l) gauge field in (2+1) dimensions 
is always confined® 2 . Because of the confinement of the 
two U(l) groups in the IGG, the remaining gauge group 
is Q = SU(2) x SU(2) x SU(2)/U(1) x U(l), where we 
identify g 1 © g 2 © g 3 with gie' 8 ^ © g 2 e i ^- 0 '> r3 © 
Particularly (—1) © 1 © 1 ~ 1 © (—1) © 1 ~ 1 © 1 © (—1) 
is the generator of the Z 2 IGG. For clarity we denote the 
element of the gauge group Q as [ 31 , 32 , 33 ]- 

We have constructed all the other 44 kinds of Z 2 spin 
liquid states with the aforementioned gauge group Q as 
showed in Table [V] The detailed dynamical properties of 
these exotic spin liquid states require further investiga¬ 
tion in the future. 


V. CONCLUSION 

In this paper, we first studied the symmetry fractional- 
ization in gapped Z 2 spin liquids on the triangular lattice, 
which is described mathematically by the projective rep¬ 
resentations of translation group, point group and time 
reversal symmetry groups. There are six independent 
gauge invariant fractionalized Z 2 quantum numbers to 
label in total 64 kinds of spin liquid states, which can 
be both illustrated using group cohomology and detailed 
analysis for triangular lattice’s symmetry group. Then 
we explicitly found that there are 20 spin liquids can be 
constructed within Abrikosov fermion formulation, while 
8 kinds can be constructed by Schwinger boson methods. 
Simultaneously, we derive vison’s fractionalized quantum 
numbers from another independent way by which vison 
in Z 2 spin liquids can be mapped to a fully frustrated 
Ising model on the dual lattice. Then we mapped all 
the bosonic states to fermionic states with the help of 






TABLE V. Exotic Z 2 spin liquid states beyond conventional 
Schwinger boson and Abrikosov fermion constructions on tri¬ 
angular lattice. Listed states No. Cl - No. C22 com¬ 
bined translational fractionalized quantum number W 12 = ±1 
give birth to all 44 states. Note that in the ip represen¬ 


tation, for example the Cl state, g M can be formally writ¬ 
ten as g M = (B^Z^gp = ir 1 © ir 2 © r° and correspondingly 


= 

sign 

m. 

i(Sm 


sign 

(nj 

uj 

)• 




No. 


10 & 

OJl 

(Ofxt 

(Oat 


9n 



S-T 


Cl 

+ 

+ 

+ 

+ 

+ 

[ir 1 

ir 2 

r U 1 

[ir 1 

ir 2 

rl 

C2 

+ 

+ 

+ 

— 

+ 

fir 1 

ir 3 

r U ] 

[ir 3 

ir 2 

r°] 

C3 

+ 

+ 

+ 

+ 

— 

pr 3 

ir 2 

r U ] 

[ir 1 

ir 3 

r°] 

C4 

+ 

— 

— 

+ 

+ 

pr 1 

ir 2 

rl 

[ir 1 , 

ir 2 , 

ir 1 ] 

C5 

+ 

— 

— 

+ 

— 

pr 1 

ir 2 

r°] 

[ir 1 , 

ir 2 , 

ir 2 ] 

C6 

— 

+ 

— 

+ 

+ 

pr 1 . 

ir 2 , 

ir 1 ] 

[ir 1 

ir 2 

r°] 

C7 

— 

+ 

— 

— 

+ 

[U 1 , 

ir 2 , 

ir 2 ] 

[ir 1 

ir 2 

r°] 

C8 

— 

— 

+ 

+ 

— 

pr 1 . 

ir 2 , 

ir 1 ] 

[ir 3 , 

ir 2 , 

ir 2 ] 

C9 

— 

— 

+ 

— 

+ 

[ir 3 , 

ir 2 , 

ir 2 ] 

[ir 1 , 

ir 2 , 

ir 1 ] 

CIO 

— 

— 

— 

+ 

+ 

pr 1 . 

ir 2 , 

ir 1 ] 

[ir 1 , 

ir 2 , 

ir 3 ] 

Cll 

+ 

+ 

— 

+ 

+ 

pr 1 

ir 2 

r°] 

[ir 3 

ir 2 

r°] 

C12 

+ 

+ 

— 

+ 

— 

pr 1 

ir 2 

rl 

[ir 1 

ir 1 

r°] 

C13 

+ 

+ 

— 

— 

+ 

pr 1 

ir 1 

r°] 

[ir 1 

ir 2 

r°] 

C14 

+ 

+ 

— 

— 

— 

pr 1 

ir 1 

r°] 

[ir 1 

ir 3 

r°I 

C15 

+ 

— 

+ 

+ 

+ 

pr 3 

ir 2 

r°] 

[ir 1 , 

ir 2 , 

ir 1 ] 

C16 

+ 

— 

+ 

+ 

— 

pr 3 

ir 2 

r U ] 

[ir 1 , 

ir 2 

ir 2 ] 

C17 

+ 

— 

+ 

— 

+ 

pr 1 

ir 1 

r U ] 

[ir 1 , 

ir 2 , 

ir 1 ] 

C18 

+ 

— 

+ 

— 

— 

pr 1 

ir 1 

rl 

[ir 1 , 

ir 2 , 

ir 2 ] 

C19 

— 

+ 

+ 

+ 

+ 

[ir 1 , 

ir 2 , 

ir 1 ] 

[ir 3 

ir 2 

r°] 

C20 

— 

+ 

+ 

+ 

— 

[ir 1 . 

ir 2 , 

ir 1 ] 

[ir 1 

ir 1 

r°] 

C21 

— 

+ 

+ 

— 

+ 

[ir 1 , 

ir 2 , 

ir 2 ] 

[ir 3 

ir 2 

T °] 

C22 

- 

+ 

+ 

- 

- 

[ir 1 , 

ir 2 , 

ir 2 ] 

[ir 1 

ir 1 

rl 


anyon fusion rules discussed previously 2 ® to give a uni¬ 
fied perspective towards parton constructions of gapped 
spin liquid states. 

Then we have carried out preliminary VMC investiga¬ 


tion on these states. We found that for J 2 /J 1 — 0.1 the 
spin order will be killed. However, in contrast to claimed 
discovery of gapped spin liquids by DMRG recently, we 
found that the Dirac spin liquid is stable for the J\ — J 2 
model on triangular lattice, and the spinon gap does not 
opened. 

In the process of constructing spin liquid states with 
slave particle approaches we found that the essential con¬ 
straint comes from the projective point gauge group, of 
which fermionic construction is SU(2) while bosonic con¬ 
struction is U(l). Motivated by this viewpoint, we pro¬ 
pose a method to construct all the other 44 kinds of spin 
liquid states beyond conventional Schwinger boson and 
Abrikosov fermion formulations by enlarging the point 
gauge group to SU(2)xSU(2)xSU(2)/U(l)xU(l). We 
also discuss one possiable realization in detail by con¬ 
structing three orbitals enbedded six spin-1/2 fermionic 
spinons to form the total spin operator within each unit 
cell. The dynamical properties of these exotic spin liquid 
states are left to future investigation. 
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